Short-range correlation in the uniform electron gas: Extended Overhauser model 
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We use the two-electron wavefunctions (geminals) and the simple screened Coulomb potential 
proposed by Overhauser [Can. J. Phys. 73, 683 (1995)] to compute the pair-distribution function 
g(r) for a uniform electron gas, finding the exact g(0) for this model and extending the results 
from g(0) to g(r). We find that the short-range (r < r s ) part of this g(r) is in excellent agreement 
with Quantum Monte Carlo simulations for a wide range of electron densities. We are thus able to 
estimate the value of the second-order (r 2 ) coefficient of the small interelectronic-distance expansion 
of the pair-distribution function. The coefficients of the small-r expansion of the spin-resolved 
So-cr'W have density or r s dependences which we parametrize in a way that makes it easy to find 
their coupling-constant averages. Their spin-polarization or £ dependences are estimated from a 
proposed spin scaling relation. 



I. INTRODUCTION 

The pair-distribution function g(r) for the uniform 
electron gas is related to the probability of finding 
a pair of electrons at a distance r from each other. 
(The expected number of electrons in a shell of volume 
4irr 2 dr when another electron is at the origin is given 
by ng(r)4:TTr 2 dr, where n = N/V is the electron den- 
sity.) In recent jy^ajcs—pajjch attention has been devoted 
to this quantityBauu'EMQ mainly because of its relevance 
in densitjznfunctional theora: g(r) is the basis of many 
non-local&B and semi-localBl exchange-correlation energy 
density functionals. A good g(r) for the uniform electron 
gas is also the starting point for building up the system- 
averaged exchange-correlation|-hole of a many-electron 
system of non-uniform densityO 

After oscillations are averaged out, the J|6«ig-range 
part (r — ► oo) of g(r) is exactly describedHO by the 
random-phase approximation (RPA) . At intermediate in- 
terelectronic distances, 0.5 < r/r s < 5 (where r s is 
the density parameter, r s = (47rn/3) -1 ' 3 ), very reli- 
able information is a.vailahle from Quantum Monte Carlo 
(QMC) simulations.BEJlalia Little is known about the 
quantitative behavior f£.g(r) for r/r s < 0.5, except in 



the high-density limit ,130 and except some estimates of 
g(0), the value of the pak^dJsbjibLLtion function at zero- 
interelectronic distanceot§EZE30'El (As the part of g(r) 
that is most transferable from uniform to non-uniform 
densities, <?(0j)-|plays a special role in spin-density func- 
tional theory£j) 

In particular, OverhauserEEl presented a simple model 
for estimating the value of g(0). The model takes into 
account two-particle correlations by means of a simple 
screened Coulomb potential with no empirical parameter. 
In his work Overhauser finds an approximate solution 
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for the dependence of g(0) on the electron density. His 
result looks realistic and performs surprisingly well in the 
known high-density limit. 

Motivated by this result, in the present work we solve 
the Overhauser model exactly, i.e., we compute the corre- 
sponding pair-distribution function for the uniform elec- 
tron gas. Comparison of our results with recent QMC 
simulations^ shows very good agreement in the range 
0-5 ;$ r/r s < 1. We thus have a strong indication that 
the Overhauser model potential gives accurate quantita- 
tive results also in the "unknown" shortest-range region 
< r/r s < 0.5. In this way we are able to present a quan- 
titative, reliable estimate for the r s dependence of the 
r 2 coefficient of the small-r expansion of g(r). This co- 
efficient is important for new energy density functionals 
which include the gradient correction to the r 2 coefficient 
of the exchange-correlation hole, while its spin resolution 
(i.e., its IT and tl contributions, also available in the 
present treatment) is of interest for functionals based on 
the Fermi hole curvature.Eil Other possible applications 
are discussed at the end of this article. 

The Overhauser model is fully quantum mechanical, 
and incorporates the effect of exchange on the pair- 
distribution function in the zero-tempexature electron 
gas. An interesting alternative approachEj mimics these 
effects by using a classical pair-distribution function with 
a non-zero effective temperature. 



II. UNIFORM ELECTRON GAS FROM 
TWO-PARTICLE WAVEFUNCTIONS 

Following Overhauser ,E£l we will construct the pair dis- 
tribution function from two-electron wavefunctions. The 
rigorous analogs of these wavefunctions are perhaps the 
"natural geiauataJs" that diagonalize the two-electron den- 
sity matrix;E3c3 the diagonal of this matrix is the pair 
density. For a generalization of the HartreerEock approx- 
imation from orbitals to geminals, see Ref.Ea. 

If we select a pair of electrons at random in the spin- 



unpolarized uniform gas, there is one chance in four that 
they will be in the singlet state, jj, — H, and three 
chances in four that they will be in one of the triplet 
states, 11, 11, 11 + |t- In the case of no electron-electron 
interaction, the corresponding two-electron spatial wave- 
functions will be 



*(r,R) = ^e 



iKR /^ik-r _i_ „-ik-r 



Eq. (^) immediately gives the exact result for a nonin- 
teracting gas, i.e., g-\i(r) = 1 for each r. To obtain the 
noninteracting g^(r) from Eq. (fil), we need to average 
over k. In the noninteracting electron gas, the probabil- 
ity distribution p{k) for k = ^|k 2 — ki| can be obtained 
geometrically by considering two three-dimensional vec- 
tors ki and k2 with < \ku2) I < kp, where 



e l r ± e * r ) , (1) Fermi wavevector. The probability p( k) is the 



where "+"is for the singlet state and "— " is for the triplet 
state, and 

R = |(ri +r 2 ), r = r 2 - ri, 

K = kx+k 2 , k=|(k 2 -ki). 

We can expand the plane waves into spherical harmonics 

oo 

e i^ = j2(2i + l)i i P i ( C os6)Je(kr), (2) 

£=0 

where Pe are Legendre polynomials and je are spherical 
Bessel functions. Then we will have 

oo 

* sing iet(r,R) = \/2V K - R J2 (2l+l)i e P e (co S 0)j e (kr) 

£ = 

even I 
oo 

*triplrt(r,R) - V2e lK - R ^(2£+l)z^P,(cos0)^(fcr). 



e=i 
oddf 



We can define spin-resolved pair-distribution functions 
for the uniform electron gas, <7i"|-(r) and g-\\_(r), corre- 
sponding to parallel- and antiparallel-spin interactions, 
and such that for the unpolarized gas 



= 2-(.9TT+5Tl)- 



(3) 



They can be related to the formulas just derived by sim- 
ple considerations: (i) if we select at random a pair of 
electrons in the uniform gas, there is a probability p(k) 
that they have relative momentum k = |k 2 — ki|/2, (ii) 
the probability for the singlet state is \ and the proba- 
bility for the triplet state is | , (iii) ^ of the triplet state 
contributes to the antiparallel-spin correlations and | of 
it to the parallel-spin correlations. So, we simply have 

9n(r) = |(l* S in g iet(r)| 2 ) + i(|* tl ,p lo t(r)| 2 ) (4) 
9n(r) = <l*tri P ict(r)| 2 ), (5) 

where () denotes average over p{k) and over the solid 
angle. Performing the spherical average over the solid 
angle, we obtain: 



g n (r) = 5>£+l)& 2 (fcr)) (6) 

OO 

g n (r) = 2£(2*+l)tff(fcr)>. (7) 
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with k ranging from to kp. (p(k) vanishes at the end- 
points of its domain, maximizes around the middle, and 
integrates to 1.) Then we find analytically that Eq. (Q) 
gives the correct result for the non-interacting uniform 
gas. Numerically, in the range < r/r s < 3 and with a 
truncation of the infinite sum over £ at ^ m ax = 7, Eq. (M) 
reproduces the known exchange-only g^ within an accu- 
racy of 10~ 6 . 

With these simple concepts in mind, we can now pro- 
ceed to compute an interacting pair-distribution func- 
tion by introducing a suitable electron-electron poten- 
tial which describes the interactions in a uniform elec- 
tron gas. Then we just have to replace the spherical 
Bessel functions ji in Eqs. (Q) and (0) with the func- 
tions Ri, solutions of the radial Schrodinger equation 
with the chosen two-body potential. Unless this poten- 
tial is very sophisticated, such a treatment will fail to de- 
scribe long-range correlations, which are mainly governed 
by collective modes, and will fail to satisy the particle- 
conservation sum rule on g aa ,(r) (Eq. (47) of Refill). 



III. SOLUTION USING THE OVERHAUSER 
SCREENED COULOMB POTENTIAL 

Overhauserta proposed a simple and reasonable model 
for the screened Coulomb repulsion V{r) in the uni- 
form electron gas. In the standard uniform-electron-gas 
model, a rigid positively-charged background maintains 
electrical neutrality. Thus Overhauser took the sphere of 
volume n _1 as the boundary within which the screening 
charge density is ne and outside of which it is zero. This 
is equivalent to assuming that the probability of finding 
three electrons in a sphere of radius r s is exactly zero, 
an assumption which is nearly true. In fact, numerical 
estimates of this probability for rap interacting electron 
gas show that it is indeed small. E3 (At r s = 5 the ratio 
between the probabilities of finding three and two elec- 
trons in the same sphere of radius r s is about 1/11; for 
larger r s this ratio is lower, and for smaller r s it is higher, 
being about 1/7 at r s = 0.) 

Thus, for interelectronic distances r < r s we expect 
the Overhauser potential to be close to the true potential 
felt by an electron moving in a uniform electron gas when 
another electron is fixed at the origin. In the region r > 
r s the potential is set to zero, and so is not expected to 
be reliable. 
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FIG. 1: Pair-distribution functions g(r\4ar the uniform elec- 
tron gas computed with the Overhauserl— I potential for differ- 
ent values of the electron-density parameter *ti Also shown 
are the Quantum Monte Carlo data from RefJ— 1. 



We also expect to have results that become more ac- 
curate as the density decreases, since the probability of 
having three electrons in the same sphere of radius r s 
becomes lower and lower. Finally, at high and interme- 
diate densities our results will be much closer to the true 
g(r) for antiparallel-spin correlations than for parallel- 
spin ones. When two electrons of opposite spins are in 
the same sphere of radius r s , a third electron is excluded 
from the sphere because of both the Pauli principle and 
the Coulomb repulsion. For a pair of parallel-spin elec- 
trons, only the Coulomb repulsion prevents a third elec- 
tron of opposite spin from entering the sphere of radius 
r s , a mechanism which becomes less efficient as the den- 
sity (and thus the kinetic energy) increases. 

In atomic units, the Ovcrhauser potential is given by 



V(r) 

V(r) 



r s \ r 







(r < r s ) 

(r > r s ). (9) 



Defining s = r/r s , q — kr s and Ut(s) = qsRi(s), the 
corresponding radial Schrodinger equation reads 



efficients c n , 

(I) _ 2 C n-i 



-re— 2 



+ r s c. 



M 



•s^-re-l 



c w -1 



(n + £){n-£- 1) 
4P =0 (n < I). 



(n>£ + l) 

(12) 



For s > 1 the solution is qs times the spherical Bessel 
function jt with a phase shift zt ■ It can be usefully writ- 
ten as 

u e (s) = qs[f e (qs)sin(qs- z e ) + (-l) l+1 f-e-i(qs) x 
cos(qs - z £ )] (s > 1), (13) 

where fj(x) are given by 



fo(x) 



h(x) 



fi-i(?) + /i+i(s) = (2j + 1) x-'fjix). 



(14) 



Equation (|13|) is properly normalized, since the noninter- 
acting case of Eqs. (||) and (0) is exactly recovered when 
Z£ = 0. By matching ui(s) and its first derivative at 
s = 1, we can find at and the phase shift zt. They are 
given by 



«7 



tan(g 



g 2 [JK g )/_ £ (g) + /i-i(g)/-i-i(g)f 

S e (q,r s )2+Ci(q,r s )2 



n / ^t s d<i> r s) 

i) - I -1 J 777 SI 

Ce(q,r s ) 



where 



Se = qf-t(q)Bt + f-t-x{q)At 
Ce = h{q)Ai - qft-i(q)Bt 

OO 

At = At(q,r s )= £ 4P(n + l) 



n=l+l 



B e = B e (q,r s ) 



n=l+l 



w 



(15) 
(16) 

(17) 
(18) 

(19) 
(20) 
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e(e+i) 



i(i+i) 



q 2 -r s 



u e = (s > 1). 



u e = (s < 1) 

(10) 



With 8 (l max = 7) p artial waves and the infinite sum over 
n in Eqs. (|ll|), ( |19| ) and (20) truncated at n max ~ 20, we 
reach accurate convergence for < r/r s < 2. 

We thus obtain smooth radial functions Ri{s^r Sl q) to 
replace the functions jt in Eqs. (0) and (R). Then we 
have to average over q to obtain g a<7 / (s, r s ). We used the 
probability distribution p(q) corresponding to an unpo- 
larized noninteracting uniform electron gas, i.e., we took 



Eq. (|) With k F = f(!% )V3 



In order to find ut(s), we proceed as follows. For < s < 
1, we expand ug in a power series: 

00 
ut(s)=at Y, c " )s " (0<s<l). (11) 

n=l+l 

Inserting this expansion into the corresponding radial 
equation, we find a recursion relation between the co- 
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P(l) = f 9 



2-(^) 1/ %+^? 3 



(21) 



with q ranging from to (97T/4) 1 ' 3 . The g-dependence 
of our results is rather weak, so we expect to have no 
significant change if we use an interacting momentum 
distribution instead of the one of Eq. ( pl| ) . 

The results for the total g(r) are shown in FisJl], to- 
gether with the newest QMC results from Ref.fLll We 



find g(r) in accurate agreement with the QMC data for 
0-5 ^ TJr a < 1 for a wide range of electron densities. In 
the shortest-range region, r/r s < 0.5, the QMC data are 
known to suffer large errors and are not so reliable. (In 
this region there is in fad_a significant discrepancy be- 
tween the data from Ref.t£l and those from RefEj.) We 
believe that for r/r s < 0.5 the present treatment pro- 
vides results much closer to the true g(r). As said, for 
r/r s > 1 the results obtained with the Overhauser po- 
tential are not reliable. 



IV. SHORT-RANGE COEFFICIENTS OF THE 
PAIR-DISTRIBUTION FUNCTION 

The small- s expansion of the spin-resolved g^ and g^ 
(where s — r/r s ), 

g n (s,r s ) = al l (r s ) + al l (r a )s + al l (r 3 )s 2 + 0($2) 

ffTT ( S ,r s ) = 4 T (r s ) S 2 + 4 T (r s ) S 3 + 0( S 4 ), (23) 

obtained by solving the Overhauser equation, has coeffi- 
cients 
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,n 



,TJ 



,TT 
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f s al l 
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(«§) 



= 6 
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' S „TT 



(24) 
(25) 
(26) 
(27) 
(28) 



where the at are given by Eq. ( |15|) and () means average 
over p(q) of Eq. (pi]). -Equations (E5T) and ( p8| ) are the 
known cusp conditions ,tHo due to the dominance for r — > 
of the term 1/r in the Schrodinger equation. 

Equation (|2J) gives the value of the pair-distribution 



function at zero-interelectronic distance, g(0) 



,TJ 



/2- 



In Fig. H we compare the present result with other 
estimates cL-.g(0): (i) two approximate solutions that 
Over hausertS obtained from his model potential after set- 
ting q = in thexadial Schrodinger equation for s < 1; 
(ii) the YasuharaLD ekctron-electron ladder evaluation; 
(iii) the Perdew-Wanga (PW92) formula, and (iv) the re- 
sult extrapolated from QMC datallS by Gori-Giorgi, Sac- 
chetti and BacheletO (GSB). The value of g(0) has been 
multiplied by the density parameter r s in order to am- 
plify the differences between the curves, which are mainly 
located at low densities. We see that the first Over- 
hauser formula, obtained from a first iteration, is quite 
far from the "exact" solution of his model presented in 
this work. His second formula, derived by a second itera- 
tion, is much closer to our solution (within 5% at r s = 2, 
50% at r s = 10). The Yasuhara electron-electron ladder 
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FIG. 2: Pair-distribution function at zero-interelectronic dis- 
tance multiplied by the density parameter r s . The present 
result (exact solution of the Ovenhpuser model) is compared 
with the two ©yerhauser formulascj (1 and 2 iterations), with 
the YasuharaHl electron-electron ladder evaluation, with the 
Perdew and Wangp (PIW92) formula, and with the extrapo- 
lation from QMC dataE3 obtained by Gori-Giorgi, Sacchetti 
and Bacheletfi (GSB). 



evaluation is the closest to our result for r s > 5 (e.g., 
the difference is only 0.6% at r s — 10). This is not sur- 
prising since both calculations are more reliable at lower 
densities. The PW92 formula is mainly a Pade approxi- 
mation of the Yasuhara values; it is not very accurate at 
densities r s > 4. where the electron-electron ladder g(0) 
goes to zero asEj rj e~ A *^ F °, while the PW92 Pade ap- 
proximant goes as rj 3 . However, the PW92 formula re- 
produces quite accurately the Yasuhara values at metallic 
densities, and has the advantage of being very simple and 
handy for many purposes. The extrapolated GSB value 
is, for r s > 0.7, much higher than the curve obtained in 
the present work, as expected from the results of Fig. [j], 
in which our g(r) curves always lie below the QMC data 
for r/r s < 0.3. 

The high-density limit of our % (r s ) is 



al 1 ^ ->0) = 1 -0.684 r s + 0(r 2 s ), (29) 

which differs from the exact one,E3E3 1 — 0.7317 r s , by 
only 6.5% in the coefficient of r s and is slightlx bet- 
ter than the— high-density limit of the YasuharaliU a , 
1 — 0.663r s .Ej We can conclude that the value of g(0) 
obtained by solving the Overhauser model is very reli- 
able: it almost recovers the exact high-density limit, and 
agrees in the low-density limit with the electron-electron 
ladder evaluation of Yasuhara. In the latter limit, where 



4 l 



and 



11 



the electrons become strictly correlatedj23 a 
tend to zero. 

Our results for g(0) can be usefully parametrized with 
the following form 
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(r.) = (1 - Br 



+ Crl + Drl 



Er 4 s )e 



-dr^ 



(30) 
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FIG. 3: Coefficients of (r/r s ) 2 in the small-r expansion of the 
pair-distribution function. The parallel- and antiparallel-spin 
contributions are separately shown. The total coefficient for 
the correlation part, i.e., a| = ^(aj +a4 ) — ^j(^p) 2 , is also 
reported. Upper panel: solution of the Overhauser model; 
lower panel: interpolation formulas including the exact high- 
density limits. 



We constrain our fit to reproduce the exact high-density 
limit, in order to be as close as possible to the true g(0). 
Thus we fix B = 0.7317 — d. A best fit of the remaining 
4 free parameters to our result for r s < 10 gives C = 
0.08193, D = -0.01277, E = 0.001859, d = 0.7524. Then 
B = -0.0207. On the scale of Fig. §, the fit error is 
invisible. 



As far as we know, the coefficients a 2 



and a 2 , given 



by Eqs. (M) and Q27j), are presented in this work for 
the first time. The good performance of g(0) obtained 
with the Overhauser potential, together with the agree- 
ment of g(r) with the QMC data for 0.5 < r/r s < 1, 
suggest that our calculation of these coefficients is quan- 
titatively reliable. In the upper panel of Fig. 0, we report 
a\ and a 2 as a function of r s . While the parallel-spin 
part shows a predictable dependence on r s , decreasing 
monotonically from the noninteracting value \{^f-) 2 ^ z to 
zero as r s increases, the || part presents at first glance 
a non-intuitive dependence on the density. However, if 



we consider the known high- and low-density behavior 
of a 2 , we see that we should expect a function which is 

qualitatively similar to the one we obtained. In fact, in 

1 1 
the r s — v limit a 2 is known to approach zero from neg- 
ative values ,13 while in the r s — > oo limit we expect a 2 
to approach zero asymptotically from positive values, in 
order to fulfill the positivity constraint on g]i(r). Thus, 



,TJ 



must have a minimum and a maximum, and must 



cross zero in between, as does our result. 

The high-density limit of our s 2 coefficients is 



(r s ^0) = -0.34 r s + 0(r 2 s ) 
(r s ^0) = i(f) 2/3 - 0.192 r s 



0(r 2 



(31) 
(32) 



Thcp.ll value is in excellent agreement with the exact 
one,Q — 0.3356 rg, while the |T coefficient of r s in tl 



11 

-0.422; see also Ref.EJ. As said, in the high-density limit 



r.h.s. of Eq. ( |32| ) ia-,54% smaller than the exact one,L 



is- 
ifl 



the assumption of zero probability of having a third elec- 
tron in the same sphere of radius r s can be nearly true 
only for antiparallel spin electrons, since in this case the 
Pauli principle keeps the third electron away. So, our a 2 
should be very close to the true one, while the ff short- 
range correlations are underestimated at higher densities 
in the present treatment (i.e., the true a 2 deviates more 
from the noninteracting value -g(^j-) 2 , so it is lower 
than our result). However, a 2 has a very simple depen- 
dence on r s , so, as we shall see, it is very easy to write 
down an interpolation formula between the high- and the 
low-density limits. 

In the upper panel of Fig. H we also show a 2 , the to- 
tal (r/r s ) 2 coefficient of the correlation part of g(r), i.e., 
\{a 2 + a 2 ) minus the noninteracting value jn(^-) 2 ^ 3 - 
We find that a 2 is positive in the density range 1.94 < 
r s < 5.88, and negative elsewhere. However, we believe 
that the positive region is a consequence of underestimat- 
ing the |T correlations. 

Our results are very well reproduced by the simple for- 
mulas 



,n 



7 TT 

'2 



(-0 a r, + laTi + S a ri + e a r 4 s ) e 

2 / 3 /i a „ ., „2\ „-d„r 



d a T s 



„. : l(9ny<* {1 _p prs + %r 2 )e 



(33) 

(34) 



where [3 a = 0.32, 7a = 0.4069, 5 a = -0.04455, e a = 
0.003064, d a = 0.4235, (i p = 0.01624, 7p = 0.00264, d p = 
0.2456. 

However, as said, the |T correlations are underesti- 
mated. We can use the formulas of Ecp. ([53|) and (J3J) to 
recover the exact high-density limits J3 by simply chang- 
ing the parameters a and d p . The new values are 
(3 a = 0.3356, dp = 0.5566. In this way, the || coefficient 
is basically unchanged, while the "["f is lowered because 
of the exact slope at r s = 0, and should be much closer 
to the true one. The new a 2 coefficients are shown in the 
lower panel of Fig. 0, together with the total coefficient 
for the correlation part. We see that now a 2 is always 
negative. 



ZiescheEJ has proposed 1 — 5(^ L ) _2 / 3 a2 l as a measure 
of correlation strength, varying monotonically from for 

r a — > to 1 for r s — > oo. 



4 r 



Here we shall motivate and use an approximate spin scal- 
ing relation for the short-range part of g: 



<w(r,r s ,C) = W(r,rr (C),C = 0), 



(41) 



V. COUPLING-CONSTANT AVERAGE 

The coupling-constant averageEil of the pair- 
distribution function for the uniform electron gas 
is given by 



g(s,r s ) = i / g{s,r' s )dr' s 



(35) 



where s is kept fixed. The function g is important be- 
cause it can account for the kinetic energy of correlation. 
Its short-range coefficients can be computed by integrat- 
ing the simple formulas of Eqs. (p0|), (|35J), and (p4|): 
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^ [(-A + Br s + Cr 2 s + Dr\ + Erf) e~ dr ° 
+A] (36) 

i [(-A! + 5 ir , + U x r 2 + ^r 3 + E ir J 
+iV s 5 )e-^+:4 1 ] (37) 

£[(-4. + S a r s + C a r 2 s + D a r 3 s + E a rf) e~ d ^ 
+A a ] (38) 



\ (x) 2/3 tJ-a* + s ^ + c ^) e ~ dprs 

+A P ], 



(39) 



where for aj' 1 we obtain A = 1.696, B = -0.2763, C = 
-0.09359, I? = 0.003837, E = -0.002471, and io±_a\ l 
we obtain Ai = 3.356, Bi = -2.525, d = -0.45, Di = 
-0.106, Si = 0.000553, F x = -0.00247. Equation @ 
has been obtained by integrating r s aj (r s ), according to 
the cusp condition. In this way, we are able to give short- 
range coefficients which exactly satisfy thejnodified cusp 
condition for g(r). (In the Perdew-Wangtl model for g, 
the cusp condition is accurately but not exactly satisfied.) 
For the s 2 coefficients, using the values corresponding 
to the exact high-density limit formulas, we have: A a — 
5.9313, B a = -2.5119, C a = -0.6997, D a = 0.03686, 



E a 



-0.007235, d a 
-0.00474, dr, ■■ 



= 0.4235, A p = 1.775, 
0.5566. 



B p = 0.0121, 



where rJ T = r s (1 + C) ^ 3 , r ^ — r s (1 — C) ^ 3 , an d 
rJ^2r s /[(l + C) 1/3 + (l-C) 1/3 ]- 

In the present treatment, information about the spin 
polarization of the system enters when we average over 
the probability p(q). For a spin-unpolarized (£ = 0) gas, 
the correct p(q) is the one of Eq. (j2~i]). For a fully po- 
larized (£ = 1) gas, we can similarly obtain p(q) from 
Eq. (ph by using the proper relation between kp and 
r g , i.e., fc^T 1 = ±{^-) 1/3 = 2 1 / 3 fc£ =0 . For a partially- 
polarized gas (0 < C < 1); we have to distinguish be- 
tween the parallel- and the antiparallel-spin cases. In the 
parallel-spin case, we can still obtain p(q) from Eq. (j^) 
by using fc FT = (1 + C) 1/3 4 = ° 

k Fi = (1 - ()^kp° 

that we can obtain pi (q) and pi (q) by rescaling p(q) of 

Eq. (Ell) as follows: 



for "ft interactions, and 
(Y^kp for J.J. interactions. This means 



plHq) = (l + (T 1/3 p[q(l + <;r 1/3 
pi l (q) = (l-()- 1/3 p\q(l-()- 1/3 



(42) 
(43) 



with q ranging from to [^(1 + C)] 1 ^ 3 f° r the TT case 
and from to [^j- (1 — C)] 1 ^ 3 f° r t ne II case. In the case of 
a pair of antiparallel-spin electrons, we have to compute 
geometrically the probability distribution for a variable 



fc = i|k 2 -k x | 
This gives 



when < |k a | < k F ± and < |k 2 | < k 



Pl l (k) 
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(forO < k< fc_) 






fc Fi) jL 1 IP i k F-j+ k 'Fl) 7„2 
5 K -f 1Z— p— p K 



" p p ' 

"Ft Fl 



(for fc_ 



"FT r 

< k< k. 



c FT . 



1 ( fc FJ+ fc Fj) 7 3 



(44) 



where fc± = i^Ft i ^-FJ,)- Our goal is to write down 
a scaling relation which allows us to use the results ob- 
tained for the C = gas. To do this, first we have to write 
down an approximate pi (q) obtained by scaling p(q). A 
reasonable choice, having correct limits and normaliza- 
tion, is 



VI. SPIN-POLARIZED GAS 

The spin-polarized electron gas is characterized by 
the parameter £ = (rij — fii)/n, where nj and n± are 
the densities of spin-up and -down electrons. The pair- 
distribution function averaged over "["f, J, J., and || pairs 
is then 

g(r,r s ,() = (^J 9n( r > r s>() + y^f) 9ii(r,r s ,C) 
+^9n(r,r s ,0. (40) 



,n 



p?(q) 



(l + f )l/3 + (1 _ C )l/3 



P 



_2q_ 



(i+C) 1 / 3 +(i-C) 1 '' 3 



(45) 



'97rNl/3 [(l+C) 1/3 + (l-C) 1/3 ] 



with q varying from to (^f ) 

In order to use the results obtained for the unpolarized 
gas, we now consider that q = kr s . Scaling q according 
to Eqs. ([42]), (|43j ) and fl45|) can thus be approximately 
equivalent to rescaling r s . This is exact if the Overhauser 
radial Schrodinger equation is changed in the same way 
when we replace q with [iq and when we replace r s with 
r s j [i. In the latter case, one also has to consider that 
the variable s = r/r s is changed into s' = ^r/r s = /i,s. 



One can easily check that this scaling property is satisfied 
by the equation with the potential set to zero (s > 1), 
which means that it is exact at the exchange-only level. 
For the case of nonzero potential (0 < s < 1), the scaling 
condition is satisfied only by the kinetic terms and by 
the 1/s part of the potential. However, since for small s 
the 1/s part of the potential is dominant, we expect to 
obtain rather good results by applying this scaling to the 
short-range coefficients. 

When we apply our scaling relation to the total s 2 
coefficient of g, we obtain 



a 2 (r 8 ,C) =(^) 2 (l + 2/3 ^,(l + C)- 1/3 ] + 

(¥)Vo 2/ M T Mwr 1/3 ] + ^ 



<- j ) 



(i+C) 1/3 +(i-C) 1/3 

2 



,n 



2/\„ 



(l +C )l/3 + (l_ C )V3 



JT 



,(46) 



where the functions a J, (r s ) and a 2 (r s ) can be cither 
the ones of Eqs. (|3^) and ( |34| ) or the ones of Eqs. J38J ) 
and (139), depending whether one is interested in g or in 
its coupling-constant average g. Equation (Eq) is exact 
(beyond the Overhauser model) in several limits: for ex- 
change, for parallel-spm correlation in the high-density 
limit (Eq. (10) of Rcf.LJ), and for the low-density limit. 
The fact that Eq. (]4£j) is exact for parallel-spin correla- 
tions in the r s — > limit suggests that the results for the 
short-range coefficients of g(r) are not affected by the vi- 
olation of the scaling relation by the part of the potential 
which remains finite when s — > 0. 

For the || interactions, we have also the error due to 
the use of an approximate relative-momentum distribu- 
tion. We can get an idea of this error by applying our 
spin-scaling relation to g(Q), with the result 

g(r = 0,r„ C) = (1 - C 2 )g(r = 0, (1+c)1/a 2 ;- (1 _ c)1 , 3 ,C = 0). 

(47) 

Equation (^) is different from-that of Marinescu and 
Quinncj or Perdew and Wang,cl obtained by using an 
exchange-like (r s — 0) ^-dependence. The high-density 
limit of g(0, r s , Q is 



g(r = 0,r a -*0,Q 

where from Eq. (|47J) 



i(l-C 2 )[l-A(C)r s 



A(C) = 



2A(C = 0) 

(i + C) 1/3 + (i-C) 1/3 ' 



(48) 



(49) 



In Fig. S-ijftf compare our Eq. (|49| ) with the exact 
A(£)/A(p)^ilH-and with the exchange-like approximation 
of Refs.cl andEj, in which this ratio is 1. We see that our 
scaling relation is always closer to the exact C-dependence 
than is the exchange-like one. The maximum error is at 
C = 1, where our ratio is 1.587 and the exact one is 1.332. 
(By the way, we numerically jConfipn the equivalence of 
the exact expressions in Refs.D andlj for all £.) 
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FIG. 4: Spin-polarization de pen dence of g(0) at high- den- 
sity. A(C) is defined by Eq. Q4q). The scaling pimriiption of 
Eq. (Ugf) is compared to the exact ^-dependenpepcJ and to 
the exchange-like scaling adopted in Refsfl andtJ. 



VII. CONCLUSIONS AND FUTURE 
DIRECTIONS 

Our work supports and extends Overhauser's picture 
of short-range correlations in the uniform electron gas of 
density n = 3/47rr^. When two electrons are separated 
by a distance r < r s , there is very little probability of 
finding a third electron within a sphere of radius r s about 
the first (especially in the strong-interaction or large-r s 
limit) . Thus the short-range correlations arise from two- 
electron collisions with an effective potential like that of 
Eq. (|). 

We have evaluated the resulting pair distribution g(r), 
not just for r = [Fig. Hi^iid Eq. (|30|), which confirm the 
accuracy of the YasuharaliU g(0)] but for all r (Fig. p. We 
find good agreement with Quantum Monte Carlo data 
in the range 0.5 r s < r < r s , where this data is accu- 
rate, while our results are probably more reliable than 
the Monte Carlo data for r < 0.5 r s . 

We have also studied the small-r expansion of the spin- 
resolved pair-distribution functions, Eqs. ( |22| ) and (|23|). 
The r s -dependent coefficients a% a (r s ) of (r/r s ) 2 from 
the Overhauser model have been extracted |-and then 
corrected for the known high-density limitQ'B [Fig. |j, 
Eqs. ( |3^ ) and (|34|) with corrected (3 a and d p ]. Unlike 



g(0) and a 2 , our a 2 has a non-monotonic dependence 
upon r s . We have also discussed how to average the 
short-range coefficients of g(r) over coupling constant 
[Eqs. (pq)-(|39|)1, and how to scale them for relative spin 
polarization £ ^ [Eqs. (Eft) and (E7J)]. In future work, 
we plan to use these results to make a new analytic model 
for the uniform-gas g(r) at all r, for use in density func- 
tional theory. Unlike currently available models which 
fail at very small or very large r s , this new model will 
be designed to work over the whole range < r s < oo, 
while satisfying more exact constraints. 



Except very close to r = 0, highly accurate pair- 
distribution functions for the uniform electron gas can 
be found from Quantum Monte Carlo simulations, or 
from the fluctuation-dissipation theorem usinga RPA* 
like response functions with the Richardson-AshcroftE-3 
local-field factor. By supplementing that approach with 
our small-r expansions of Eqs. (pa) and (|23|), it should 
be possible to find essentially-exact spin-resolved pair- 
distribution functions and correlation energies. 

Finally, we suspect that this approach can be applied 
usefully to Ojtiiei; systems, including the two-dimensional 
electron gasSEa 
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